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Recently, Meyer and Wallach [D.A. Meyer and N.R. Wallach (2002), J. of Math. Phys., 
43, pp. 4273] proposed a measure of multi-qubit entanglement that is a function on pure 
states. We find that this function can be interpreted as a physical quantity related to the 
average purity of the constituent qubits and show how it can be observed in an efficient 
manner without the need for full quantum state tomography. A possible realization is 
described for measuring the entanglement of a chain of atomic qubits trapped in a 3D 
optical lattice. 

Keywords: Entanglement 

1 Introduction 

Characterizing many body entanglement has become an active area of research both for under- 
standing dynamics of interacting systems and identifying resources for quantum information 
processing (QIP). Many tasks in QIP involve mapping a configuration of spins prepared in a 
product state to a number of different classes of entangled states. In order to characterize the 
dynamics of entanglement as well as the performance of the mapping, it would be helpful to 
have a single parameter that quantifies the amount of multi-particle entanglement contained 
in the system at any given time step. A good measure of entanglement would be one that 
faithfully captures the nature of nonlocal quantum correlations in a quantum state and as 
such should be a mapping from states to reals that is non increasing under local operations 
and classical communication on average. Such measures are known as entanglement mono- 
tones |2| . There is no unique function that satisfies these properties, and for multi-particle 
systems the problem is compounded by the fact that different degrees of entanglement can be 
shared between different partitions of the total system. Several important results have been 
found. For example, for bipartite pure states, a family of entanglement monotones is given 
by the Schmidt coefficients of the quantum state |2]. For more than two subsystems, various 
functions have been introduced to quantify entanglement. These include the fully entangled 
fraction which measures the ability of a state to perform tasks such as teleportation and 
dense coding, polynomial invariants, which involve moments of the reduced state eigenvalues 
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0] , and the Schmidt measure 5 , which is related to the minimum number of terms in the 
product state expansion of a state. The later is an entanglement monotone, as are some 
polynomial invariants. 

2 Pure state entanglement 

A difficulty with many entanglement measures is that they do not directly correspond 
to the outcome of a set of observables on the state. It has been shown that with a small 
number of identical copies of a state, one can measure quantum entanglement between two 
systems without prior knowledge of the state For more than two systems, most entangle- 
ment measures require knowledge of the state itself, which involves performing quantum state 
tomography. For a state residing in a d dimensional Hilbert space, this means a minimum 
of d 2 measurements for each statistical sample must be made on the system. It would be 
advantageous, especially with regard to long spin chains, to have an easily measured quantity 
of multi-qubit entanglement. Meyer and Wallach 1 have proposed a multi-qubit entangle- 
ment measure which was derived as a formal function on many qubit states but has a simple 
physical interpretation. They introduce a function on pure states of n qubits defined as, 

A n 

Q(|V>>) = -££(|u fe >,|S fe >), (i) 

fe=l 

where \u k ) and \v k ) are vectors in C n which are non-normalized (indicated by the") and 
obtained by projecting on state \ip) with local projectors on the fcth qubit, 

IV) = \0}k®\u k ) + \l) k (g>\v k ). (2) 

The function D(\u k ), \v k }) measures a "distance" between the two vectors \u k ) and \v k ). It is 
obtained by taking the generalized cross product: 

D(\u k ),\v k )) =J2 \u^ -u^\ 2 . (3) 

The quantity Q(\ip}) is invariant under local unitaries because the effect of applying Uj is to 
multiply each term D(\u k ), \v k }) by \det(Uj)\ 2 = 1. One could also write the state \tp) in the 
Schmidt decomposition over the bipartite division of qubit k and the other qubits as: 

m = \0) k ®\i k ) + \l) k ®\y k ) 7 (4) 

where (x k \y k ) = 0, and {|0)fc, |l)fc} are related to {|0)fc, |l)fc} by a local unitary operator U k . 
The purity of the state of qubit k is therefore Trip 2 } — (x k \x k ) 2 + (y k \y k ) 2 ■ By invariance 
of the generalized cross product under local unitaries, D(\u k ), \v k )) — D(\x k ), \y k )) and this 
quantity can be related to the norm of an antisymmetric tensor M k — \x k )(y* k \ — \y k )(x* k \ 
as [7] 

D{\z k )Ay k )) = E^-lW-fiKl 2 

= \Tr[MlM k ] (5) 
= (x k \x k )(y k \y k ) 
= 1/2(1 -Tr[p 2 k }). 
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Therefore, 

n-l 

Qm=2(l-l/nJ2Tr[p 2 k }). (6) 

fc=0 

This expression elucidates the physical meaning of the multi-particle entanglement as an 
average over the entanglements of each qubit with the rest of the system. Put in this form, 
the entanglement measure can also be identified with the purity of a state relative to the Lie 
algebra of local observables for n qubits within the generalized-entanglement framework of 
Barnum et al. 

The measure Q is linearly related to the average purity of each qubit and it is easy to 
check that it satisfies the following properties: (i) < Q(\ip)) < 1. Q{\ip)) — iff \ip) is a 
product state and Q(\ip)) = 1 for some entangled states, (ii) Q(\i/j)) is invariant under local 
unitaries Uj. 

Property (i) follows from the fact that any reduced state, pj, is at most rank two and 
because its eigenvalues sum to one, 1/2 < Tr[pj] < 1. When \ip) is a product state, Tr[pj] = 
1 Vj implying Q{\^)) = 0. Conversely, Q(\ip)) — implies X^=o Tr\pf\ = n but because 
Tr[p"j] is bounded above by 1 this can only be true if Tr[p 2 ] = 1 V). Q{\ip)) — 1 only if 
Tr[p 2 ] = 1/2 Vj meaning the reduced state of every qubit is maximally mixed. Property (ii) 
is a consequence of the fact that the eigenvalues of the reduced states are invariant under 
local unitaries. 

Many pure states obtain the maximum value of Q including the n-partite GHZ state: 
l/v^dO . . .0) + |1 . . . 1)), and the "cluster state" [H] defined as the n qubit state: 

n— 1 

|*« clus) = (g)(|0)a^ +1 + |l)a). (7) 

a=0 

with the convention <r™ = 1. These two states are inequivalent under local unitaries as best 
seen by comparing the Schmidt number of each, defined as the minimum number of product 
states over all subsystems necessary to describe the state. In the former it is two, while for 
the later it is 2 L"-/ 2 J , it is relatively easy to calculate the entanglement of certain pure states. 
For instance, the n qubit PF-state ^U] defined as \W n ) — ^= l 71 "^ • • ■ 01)): where the sum 
is over the n cyclic permutations, has an entanglement value Q(\W n )) = 4(n — l)/n 2 . Meyer 
and Wallach 1 calculate the entanglement of several other states such as the eigenstates 
of the Heisenberg interaction and single logical qubit states encoded in the five qubit error 
correction code. 

A deficiency of Q as an entanglement measure is that it cannot distinguish sub-global 
entanglement. For example, the two states |$) = 1/V2(|00)i, 2 + |H)i,2) ® 1/V2(|00) 3: 4 + 
1 11>3,4> and |*) = 1/V2(|0000)i, 2 ,3,4 + |llll)i,2,3,4) both have Q values equal to 1. One could 
distinguish different states with the same Q value by comparing the Schmidt numbers over 
various bipartite divisions. For instance, the Schmidt number of |$) over the division into 
two sets of qubits {1,2} and {3,4} is 1, while for |*I>) it is 2. 

2.1 Observation of entangled states 

A benefit of the reformulation of measure Q in Eq. [fj] is that it identifies a relatively small 
number of parameters that need to be estimated in order the measure the multi-qubit en- 
tanglement. The computation of an arbitrary function on a state of n qubits requires good 
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estimates of all 2 2n parameters of the state. By the formulation of Q in terms of the average 
purity of the qubits, it is only necessary to estimate a number of parameters linear in n, in 
particular only the purity of each qubit. 

The purity of a single quantum state can be observed using an entanglement witness as 
proposed by Filip IT]. He shows that given two systems, A and B, initially in the product state 
p(E> p, one can measure a quantity proportional to Tr[p 2 ] by interacting the two systems with a 
single ancillary qubit and measuring the ancilla. The ancilla is prepared in the superposition 
l/v^QO) + |1)) and plays the role of the control qubit in a controlled SWAP (c-SWAP) gate 
acting on the two systems. When the ancilla is measured in the a x basis, the probability to 
measure population in the ± eigenstate is: p(±) = 1/2(1 ±Tr[p 2 ]). In ^JJ it is claimed that 
this constitutes a quantum non-demolition (QND) measurement of purity of the quantum 
state p because the output of the quantum process is a balanced mixture of the two quantum 
states regardless of the measurement outcome. This, in fact, is not the case because the 
systems generally become entangled with each other. The measurement process described 
above constitutes a measurement of the SWAP operator which projects the two systems into 
the symmetric or antisymmetric subspaces of the input state p ® p. Specifically, if we write 
the state p in an eigendecomposition: p = ^iP^i)^^, then a selective measurement of the 
ancilla in the ± eigenstate of cr x corresponds to the map: 

p®p^p'{±) = ——Y^PiP j {\(l>i)A\<t>j)B± |^)a|^)b)U(^|b(0j| ± A (<Pj\B((t>i\). (8) 

The reduced states of the output state p'(±) are again identical but now equal to: Tr^[p'(±)] = 
Ttb[p' '(±)] = (p± p 2 )/(2p(±)). Because p'(±) is an eigenstate of SWAP, subsequent mea- 
surements using the entanglement witness will reveal no further information about the state 
overlap Tr[p 2 ]. Further, without knowledge of the input state it is not possible to unitar- 
ily disentangle the two systems. However, a significant advantage is that the measurement 
consumes only two copies of the quantum state. This should be compared to quantum state 
tomography onaci dimensional system where d 2 — 1 measurements are needed to reconstruct 
the state and the same number of copies are consumed. In the case of qubits, a measurement 
of purity using the entanglement witness is 50% more efficient than performing state tomog- 
raphy. For both techniques, the measurement statistics on N identically prepared ensembles 
converge to the correct value of purity with an error that scales like 0(l/\/~N). 

In order to measure Q for the state of a chain of spins, it is sufficient to measure the purity 
of each spin. A schematic for implementing such a measurement in an efficient, parallel manner 
is shown in Fig.^ The setup includes a control register c and two identically prepared target 
registers t and s all of length n. It is assumed that interactions between qubits within a 
register can be suppressed during the measurement process so that interactions take place 
only between triplets of qubits {cj, tj, Sj}, where j is the site index. The c-SWAP gate can be 
implemented in parallel over the columns of qubits with only nearest neighbor a z ® a z type 
interactions between pairs (cj, tj) and (tj, Sj) as we now show. 

In deriving a pulse sequence for c-SWAP, we first point out that the SWAP gate, which 
is generated by the Heisenberg interaction H exc = ga ■ <?, can be simulated by the Ising 
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interaction, H zz = ga z ® <j z , by conjugation with single qubit gates as: 
SWAP(i, s) = e i ^e~ i ^^® a * +a y® a y +at z® aB * ) 

This simulation minimizes the time spent during two qubit interactions |12|. To implement 
the c-SWAP gate, we replace e - ^^®" 7 * in the above sequence with e - i il 1 ><==( 1 l®< 7 * ! ® 1 = 
g-'-g- o- i «io- Z g4-g-(T a .«i(T i «io- z _ rpjjg three body interaction can be simulated with pairwise interac- 
tions as, 

The total sequence for the c-SWAP is therefore: 

c - SWAP(c, t, s) = e _i f ^e^^+^^e^f CT ^® CT *e i t CT J®^® CT ie~ i 5^ +CT ^e^(°'* +<7 ^ 

g— i-f-cr* <g><T° gj^ cr°(g><r* ®<7° g— i-J (c* +o-°)g — * J <T z®< T z gj-l± a 'z® a \® a % 

(11) 

For two qubit gates generated by pairwise Hamiltonians, i? zz , with g > 0, the total 
interaction time needed to implement the c-SWAP gate using this sequence is T c _swap = 
277r/(4g). If the sign of g is tunable to g < 0, then the interaction time is reduced to 
T c _swap = 97r / (4 1 (7 1 ) . After performing c-SWAP gates in parallel, a measurement of the 
qubits in register c in the cr x basis will yield a mean number of l's equal to: n(p(+)) = 
nQ(M)/4. 
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Fig. 1. Stacked array of qubit registers as a platform to measure entanglement properties. Two 
registers of n spins each, denoted t and s, are identically prepared in the state \tp) and interact with 
a third register of qubits c prepared in [l/v^2(|0) + A sequence of gates implements the 

controlled-SWAP gate in parallel over the columns, with the qubits labelled Cj acting as controls 
to the targets tj , Sj . After the sequence, the c register is measured in the a x basis to extract 
information about the entanglement of 

It will be important to verify that the state of a system s is a pure state. If it can be 
assumed that the two systems s and t are identically prepared, each in state p, then one can 
measure the purity of the entire state by a simple adaptation of the above procedure. If the n 
bit control register c is prepared in the initial state |</>(0)) = 1/V2(|0) + |1)) <£> |0 . . . 0), then 
a sequence of gates U = ®i._ n _i CNOT(/c — 1, k) will produce the n-partite GHZ state in c. 
Now the c-SWAP gate over the entire states s and t can be applied bit-wise over the triplets 
of qubits (cj, tj, Sj). After this the control register is disentangled by operating with U^ 1 and 
the qubit cq is measured in the a x basis to measure the purity of the state p. Using entangled 
states of the control system allows a measurement to infer the purity of any subset of qubits in 



624 An observable measure of entanglement for pure states of multi-qubit systems 

the state of s simply by preparing an entangled state of the corresponding control qubits and 
applying the c-SWAP gate over that subset. In this way, states with sub-global entanglement 
can be distinguished from globally entangled states not separable over any subsystem. After 
measuring the purity of a given subset of qubits, the constituent qubits in the two registers s 
and t will become entangled with each other. In order to disentangle them and reuse them, the 
state of the qubits in the measured set of the t register will need to be depolarized and reset to 
|0). In general, unless there exists prior knowledge of the distribution of entanglement across 
the system, this technique for distinguishing different partitions of entanglement is inefficient 
because the number of partitions of the n qubits is 2 n . 

A testbed for implementing this multi-qubit entanglement measurement is an ensemble of 
trapped neutral atoms in a three dimensional optical lattice. This system has the advantage 
of inherent parallelism for quantum processing over identically prepared registers of atomic 
qubits. There are several proposals for QIP in optical lattices with varying demands on 
addressability of the atomic qubits. Because Q is a quantity related to an average signal over 
the control register, some constraints on addressability within the system can be relaxed. If 
it can be safely assumed that state of the register is pure, then by using the three neighboring 
lattices as explained above, one does not need to resolve qubits within each register but 
only between registers c, t, and s. Therefore, the following protocol should be applicable 
both to implementations using a polarization gradient lattice along one dimension where the 
inter-atomic spacing is not directly optically resolvable and schemes where the atoms 
are trapped along that dimension by long wavelength light of uniform polarization |14l . 
Atoms can be trapped along the other two directions by two pairs of counter-propagating 
linearly polarized beams that are distinct in frequency so that they do not interfere with 
each other. Such a configuration produces planes of identical one dimensional registers with 
planar separation along the x direction. In order to measure the multi-qubit entanglement in 
a given register, it is necessary to have a mechanism to interact neighboring planes of qubits, 
and have an inter-plane separation much larger than an optical wavelength so that the planes 
can be addressed. A possible solution is to use a bi-periodic potential to trap along x as is 
suggested in 15 . There it is proposed to trap 87 Rb atoms along one dimension with two 
pairs of counter-propagating linearly polarized laser beams from the fundamental and first 
harmonic of a CO2 laser. By adjusting the relative intensities of the two colors of light, pairs 
of atoms can be made to interact along that dimension. 



Fig. 2. Cross section of 3D optical lattice shown as a platform to implement a measurement of 
entanglement. The trapping potential along x results from an incoherent sum of two light fields 
and can be tuned to join planes of trapped atoms, ordered (c, t, s), together pairwise. (a) Initial 
configuration with Vi = 0, 5 = 0. (b) Vl/Vq = 0.3, <5 = 0. (c) Vb/Vi = 0.3, <5 = vr/2. 

Consider three planes of atoms, ordered (c, t, s), localized along x at the minima of an 
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external potential V(x) — Vocos 2 (kx) with inter-planar separation Sx = w/k. The planes 
could each constitute a single quantum register, or an array of ID registers. Introducing a 
second light field, phase locked with respect to the first by a phase S, and at half the frequency 
creates a bi-periodic potential: V(x) = Vo cos 2 (kx) + V\ cos 2 (kx/2 + 5). If the second field 
is turned on adiabatically with respect to the oscillation period of the wells created by the 
first potential, then the motional states of the atoms will follow the ground state of the new 
potential. As shown in Fig. by adjusting the relative intensity and phase of the two fields, 
atoms in planes (t, s) or (c, t) can be joined pairwise at the two minima of a double well 
potential. The double well barrier height can be adjusted to allow substantial overlap in 
the motional wavefunctions of each atom in neighboring planes. If the quantum information 
is encoded in the proper ground state electronic or motional levels, a controlled ground- 
ground state collision can be engineered to generate evolution by an Ising type interaction 
|17|. Provided the fundamental trapping frequency is chosen so that Sx > Tr/k res , where 
ck res / (2tt) is the resonant atomic frequency, the 3 planes of atoms will be addressable. Because 
the potential period actually depends on the difference wave- vector of the counter-propagating 
beams, addressability can also be mandated by using beams that are not parallel, but rather 
intersect at an angle 8 so that Sx is increased by a factor of 1/ sin(0). The c-SWAP gates can 
be implemented in parallel and the values of the qubits in the control register can be read out 
at once using a spatially broad laser on a resonant cycling transition. The observed resonant 
fluorescent signal will be proportional to the average purity. 

3 Conclusions 

We have shown that a previously introduced entanglement measure on pure states of an 
arbitrary number of qubits is simply expressed as a linear function of the average purity 
of the qubits. This quantity can be observed in a straightforward way in the laboratory 
either by performing state tomography on each system qubit or by a more efficient technique 
using bitwise interactions with a second identically prepared register. The latter protocol is 
particularly well suited for systems for QIP with massive inherent parallelism such optical 
lattices. It is hoped that other, directly observable, measures will be found that more sharply 
distinguish the partitioning of shared entanglement over multi-partite systems. 
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